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Abstract 

We derive explicit semiclassical quantisation conditions for the Dirac and Pauli 
equations. We show that the spin degree of freedom yields a contribution which 
is of the same order of magnitude as the Maslov correction in Einstein-Brillouin- 
Keller quantisation. In order to obtain this result a generalisation of the notion of 
integrability for a certain skew product flow of classical translational dynamics and 
classical spin precession has to be derived. Among the examples discussed is the 
relativistic Kepler problem with Thomas precession, whose treatment sheds some 
light on the amazing success of Sommerfeld's theory of fine structure [Ann. Phys. 
(Leipzig) 51 (1916) 1-91]. 
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1 Introduction 



Semiclassical quantisation conditions provide the most direct link between the old quantum 
theory of Bohr and Sommerfeld on the one hand and wave mechanics on the other hand. 
Unlike other semiclassical tools, like trace formulae or the Van Vleck-Gutzwiller propagator 
||l|, they do not only express quantum mechanical objects in terms of classical properties 
but also employ exactly the same formulation as was used by the old quantum theory, 
namely action quantisation. 

Before the advent of quantum mechanics quantisation of a system was done by deter- 
mining action integrals within the classical theory and setting these equal to an integer 
multiple of Planck's constant h = 27ih, i.e. one required 

^ pdx = 27ihn (1-1) 

with integer n. This condition, originally put forward by Bohr ||^ in order to understand 
the hydrogen spectrum, was first understood as the quantisation rule for one degree of 
freedom. 

Around 1915 there was an ongoing discussion how this condition should be translated 
to more than one degree of freedom, see e.g. |^, ||, |^ and the introduction of [If. Epstein 
proposed to use that set of coordinates in which the problem separates if such coordinates 
exist. For each degree of freedom there would then be a condition of the form ([0| ), i.e. 

j) Pj dxj = 2'Khnj (1.2) 

with integers n^, j numbering the degrees of freedom and (pj, Xj) being a pair of canonically 
conjugate variables in that particular set of coordinates. Epstein's point of view was 
assumed by Sommerfeld who successfully applied this prescription treating many problems 
in spectroscopy [|]. 

Shortly after, Einstein pointed out that separability of the equations of motion is 
not a necessary condition for action quantisation but that merely integrability (in the sense 
of Liouville and Arnold 0, 3) is required: If there are sufficiently many integrals of motion 
with pairwise commuting Poisson brackets then the phase space foliates into invariant tori 
on which the line integral 

Xf 

pAx (1.3) 
is locally path-independent. The quantisation conditions can then be written in the form 

® pdx = 27rhnj (1.4) 

where {Ci, . . . ,Cd} denotes a basis of non-contractible loops on a given torus. This for- 
mulation has the advantage over Epstein's version of being independent of the coordinate 
system. 
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Soon after the old quantum theory had been replaced by matrix and wave mechanics 



the old quantisation conditions were rederived and modified by Wentzel ||10l, Kramers pi 
Brillouin flT 



and Jeff'reys ||T^ in a short-wavelength approximation, the so-called WKB or 
JWKB method. It was shown that depending on the character of the motion the quantum 
numbers may have to be shifted by a small number, leading e.g. to half-integer quantum 
numbers for oscillations but integer quantum numbers for rotations. Again the treatment 
was first for one-dimensional and then for separable systems. A complete derivation of 
the quantisation conditions from the Schrodinger equation that takes into account both 
the abstract integrability condition used by Einstein and the small shift of the quantum 
numbers is due to Keller |M. He proved the semiclassical quantisation conditions 



pdx = 2'Kh [nj + ^ 



(1.5) 



which are now known as Einstein-Brillouin-Keller (EBK) or torus quantisation. The num- 
ber fij G Z4 denotes the Maslov index, see e.g. ||1^, |T6|, a topological invariant of the cycle 



Cj. In one dimension it counts the number of turning points encountered along the loop. 



A good overview on theses topics is given in ||17 



The discussion so far was only for non-relativistic quantum mechanics, i.e the semi- 
classical approximation for the Schrodinger equation. In 1916 Sommerfeld applied the 
quantisation conditions (|1.2| ) also in a relativistic context 0. His aim was to find small 
corrections to the hydrogen spectrum which he expected to be due to relativistic effects. 
The success was overwhelming, the so-called Sommerfeld fine structure formula agreed 



excellently with the experimental data. More than ten years later [T^, [T^ it was found 
that the energy levels of the hydrogen atom, when calculated using the Dirac equation, the 
correct relativistic wave equation for the electron, are identical to the levels determined 
by Sommerfeld. The Dirac equation, however, does not only take into account relativistic 
effects but also the half-integer spin of the electron. The fine structure accounted for by 
the Sommerfeld formula is to a large extent due to spin-orbit coupling, an effect that was 
unknown at the time Sommerfeld did his calculations. In fact, even the property of spin 
itself was yet to be discovered. This seeming paradox has to be explained by a semiclassical 
analysis of the Dirac equation. 



Early semiclassical approaches to the Dirac equation are due to Pauli |22| and Rubinow 
and Keller ||21|. These will be discussed in section ^ The subtleties connected with such 
an approach are related to the fact that the Dirac equation is a partial differential equation 
for a spinor and not just for a scalar wave function. Therefore these problems should be 
discussed in the more general context of semiclassical (or short wavelength) approximations 



to multicomponent wave equations. It was observed by Yabana and Horiuchi ||22| that the 



occurrence of geometrical or Berry phases [^, |2J| plays an important role in this context. 
Kuratsuji and lida ||2^, |26|, using path integral methods, suggested that the symplectic 



structure of phase space should be deformed such that it includes the contribution of 



geometric phases. A review of these results is given in ^j. A general method for the 
semiclassical quantisation of multi-component wave equation was derived by Littlejohn and 
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Flynn ^7j. Their results hold whenever the principal Weyl symbol of the Hamiltonian, 
which is a matrix valued function on classical phase, has non-degenerate eigenvalues. Thus 
it does not apply to the Dirac equation, in which case these eigenvalues have a multiplicity 
of two as will be shown in section ||. In such a situation the geometrical phases that 
appear in the semiclassical expressions are not simple U(l)-phases any longer but are 
themselves matrix-valued and thus in general non-Abelian. It was shown by Emmrich and 



Weinstein ||2^ that in this case integrability of the dynamics generated by the classical 



ray Hamiltonians (the eigenvalues of the principal Weyl symbol) is no longer sufficient to 
guarantee the existence of global semiclassical wave functions and thus of semiclassical 
quantisation conditions. 

In this article we derive semiclassical quantisation conditions from the Dirac equation, 
the relativistic wave equation for particles with spin |, and from the Pauli equation, de- 
scribing particles with arbitrary spin in a non-relativistic context. We show how to resolve 
the problems mentioned in the context with the occurrence of non-Abelian Berry phases 
by developing a generalisation of the notion of integrability that not only includes the 
dynamics of the ray Hamiltonians but imposes an additional condition. In this way we 
effectively reduce the non-Abelian phases to U(l)-phases. The latter enter the semiclas- 
sical quantisation conditions by a correction which is of the same order as the Maslov 
contribution in ( |1.5|) . This correction represents the influence of the spin degree of freedom 
and can be given a clear physical interpretation in terms of classical spin precession. By 
applying the method to the relativistic Kepler problem we shed some light on the success 
of Sommerfeld's fine structure formula. A brief account of some of these results was given 
in 
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The paper is organised as follows. In section | we outline the derivation of EBK- 
quantisation for later reference, thereby emphasising the role of integrability. Section ^ 
deals with the determination of semiclassical wave functions for the Dirac equation. In 
section || we generalise the concept of integrability to the case of group extensions and the 
skew product of classical translational dynamics and classical spin precession. Based on 
this characterisation we then derive explicit semiclassical quantisation conditions for the 
Dirac equation in section |. In section ^ we show how these results translate to the Pauli 
equation with arbitrary spin. Before we treat some special examples in sections || and ^ 
(among which is Sommerfeld's theory of fine structure) we derive general formulae which 
facilitate the semiclassical quantisation of spherically symmetric systems in section |^. We 
conclude with a summary in section Some important formulae for Weyl quantisation 
are listed in appendix |^. 



2 EBK quantisation 

In this section we briefiy summarise the main steps in the derivation of the EBK quantisa- 
tion rules. We do so in order to introduce some notation and for later reference such that 
we can explicitly compare to this basic situation when treating systems with spin. This is, 
however, not intended to be a complete review of EBK quantisation and we thus refer the 
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reader seeking a comprehensive introduction to EBK quantisation to the hterature, e.g. 



14, 17 



We want to find an asymptotic solution of the stationary Schrodinger equation, 

HiIj{x) = EiP{x) , (2.1) 

where G L^(M'^) and H shall be a Weyl operator (some facts on Weyl quantisation are 
summarised in appendix ^ with symbol 

H{p,x) = Ho{p,x) + hHi{p,x) + {h'^), h^O. (2.2) 

The leading order term Hq in the semiclassical limit h is known as the principal symbol 
of H and Hi is called the sub-principal symbol. For simplicity we will choose Hi = for 
the rest of this section. For the familiar Hamiltonian 

H = -—A + Vix) (2.3) 
2m 

describing a particle of mass m moving under the influence of the external potential V{x) 
the Weyl symbol reads 

H{p, x) = Hoip, x) = ^ + V{x) (2.4) 

and we can write H = iJ(y V, x). 

For the wave function ip one tries the WKB ansatz 

i^WKB{x) = J2(f) a,{x) e-^'^'^K (2.5) 

fc>0 ^ ^ 

Inserting into ([2?l] ) yields in leading orders in h, cf. appendix 0, 

h 



[H{V^S, x)-E]ao + j[ [H{V^S, x) - E] ai 
+ [{Vj,H){V^S, x)] (V^ao) + ^ [V.(Vpi/)(V.5, x)] ao} + {h') = . 



(2.6) 



This equation is easily confirmed for the particular Hamiltonian ( |2.3| ) by direct computa- 
tion, but it also holds for arbitrary (semiclassical) Weyl operators, see appendix 0. The 
strategy of the WKB method is now to satisfy eq. (|2.(]| ) separately order by order in h. In 
leading order one finds the Hamilton-Jacobi equation 

H{V^S,x)=E (2.7) 

of classical mechanics with the (principal) symbol H{p, x) acting as the classical Hamil- 
tonian. From standard Hamilton-Jacobi theory, see e.g. ||3T|, ||, one thus concludes that 
the phase S{x) of the WKB ansatz becomes the classical action generating the dynamics 
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with Hamiltonian H{p, x): If X(t)) is a solution of Hamilton's equations of motion 

then VxS{X{t)) = P(t). A solution of the Hamilton- Jacobi equation ( p.7|) can thus be 
obtained by integration along solutions of Hamilton's equations of motion in the following 
way. Denote by y an arbitrary point in configuration space and by ^ a momentum satisfy- 
ing H{$,, y) = E. Let 0^ be the (Hamiltonian) fiow generated by the classical Hamiltonian 
H{p, x)^ i.e. 0^(^, y) =: {P(t), X{t)) denotes the point reached at time t on the trajectory 
starting at {^,y). Then we have 



Six) - Siy) 



dt' 



VxS{X{t'))X{t')dt' = / PdX (2.8) 

Jy 



where in the last expression integration is along the trajectory 0^(^, y). Finally the action 
reads 

S{x) = S{y)+ I PdX, (2.9) 
-'y 

where S{y) is the arbitrarily chosen value of S at the point y. Given a solution (|2.9| ) of the 
Hamilton- Jacobi equation the next-to-leading order equation deriving from ( |2.6| ) reduces 
to 

[(Vpi/)(V.S, x)] (V.ao) + ^ [V,(Vpif)(V,5, x)] ao = . (2.10) 

This is known as the transport equation for the leading order amplitude ao{x). Due to 
Hamilton's equations of motion the first term can now be interpreted as a time derivative 
along the the trajectory 0//(^, y) which we shall denote by ^ or simply by a dot, 



a = ^:=[{VpH){VxS,x)] (V.a). 



(2.11) 



The solution of (|2.10|) is locally given by 



ao[x) 



det 



dy_ 

dx 



(2.12) 



see e.g. 14, 17 . Together with (p.9|) one has thus found an approximate solution 



1 dy 

det 7— exp 
ox 



PdX 



) 



(2.13) 



of the Schrodinger equation for points in a neighbourhood of y that are visited along a 
a solution of Hamilton's equations of motion starting at y. 

However, one still needs to find a way to integrate the Hamilton- Jacobi equation along 
paths transversal to 0^. Furthermore, the approximation ( |2.13|) breaks down at points 
where 0//(^, y) touches a caustic and thus || becomes singular. For one degree of freedom 
these points are given by the turning points. In order to address both these problems one 
has to know more about the classical phase space structure. For arbitrary Hamiltonians one 
can in general not proceed far beyond this point as was already pointed out by Einstein 
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101 in the context of the old quantum theory. Instead one has to invoke the concept of 
integrability. 

Following Liouville we say that a Hamiltonian is integrable if there are d constants 
of motion Ai := H, A2, . . . ,Ad which are in involution, i.e. whose Poisson brackets vanish 
pairwise. A complete proof of the consequences of this definition is due to Arnold whose 
version chapter 10] we shall quote here. 

Theorem 1. (Liouville- Arnold) Suppose that we are given d functions in involution on a 
d- dimensional symplectic manifold 

{/l,-,Afc} = 0, j,A; = (2.14) 

Consider a level set of the functions Aj, 

Ma = {(p, x) I Aj{p, x) = aj, J = l,...,d}. (2.15) 

Assume that the d functions Aj are independent on Ma (i.e. the d 1-forms dAj are linearly 
independent at each point of Ma). Then 

1. Ma is a smooth manifold, invariant under the phase flow with Hamiltonian function 
H = Ai. 

2. If the manifold Ma is compact and connected, then it is diffeomorphic to the d- 
dimensional torus 

T"^ = {(dj,...,'dd) mod 2n}. (2.16) 

3. The phase flow with Hamiltonian function H determines a conditionally periodic 
motion on Ma, i.e. in angular coordinates i? = (dj, . . . ,^d) we have 

— = u), uj = uj{a). (2.17) 

dt ^ ^ ^ ' 

4. The canonical equations with Hamiltonian function H can he integrated by quadra- 
tures. 

We refrain from providing a proof of the theorem here but refer the reader to Arnold's 
book Q. Instead we remark on some aspects which are relevant for the following sections. 

In order to prove property 2 one first shows that the conditions ( |2.14| ) imply that the 
fiows 0* generated by the observables Aj, j = 1, . . . ,d, commute on Ma, i.e. 

</>5° = 4' yj,k = i,...,d. (2.18) 

this yields a transitive action of M*^ on Ma. In the following we will refer to Ma as a 
Liouville-Arnold torus. 

In general the coordinates i? and the observables A are not canonically conjugate. 
However, (locally) there exists a mapping A I such that (7, i?) form a set of canonically 
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conjugate variables. The new constants of motion I{A) are called action variables and 
the explicit construction of action and angle variables (/, t?) is the desired integration of 
of Hamilton's equations of motion by quadratures. The Hamiltonian expressed in the new 
variables becomes a function H{I) of the action variables only and thus the frequencies in 
( p.lTp are given hj uj = V/iJ(/). 

Theorem |I] can be used in order to derive the EBK quantisation conditions as follows. 
Since the flows 0* commute we can deflne the action S{x) analogously to ( p. 9]) by integra- 
tion along the flow lines of (^2, ■ ■ ■ , 4>d instead of (f)\ = 0^. To this end define the multi-time 
flow 

$* := 0*^* o . . . o , (2.19) 
where due to (|2.18| ) the ordering is unimportant. Since is a transitive action of M"^ on a 



Liouville-Arnold torus, for any a; in a small neighbourhood of y there is a unique t e M 
such that 

^\^,y) = ip,x) (2.20) 

with some momentum p. The rapidly oscillating phase of the WKB wave function is then 
given by 

S{x) = S{y)+ [ PdX (2.21) 



where integration is along $*(^, y). Since the flows 0*, j = 1, . . . , (i, commute this is not 
in conflict with the requirement of S solving the Hamilton- Jacobi equation ( p.7|) . This 
answers the question of how to integrate transversal to 0^. 



One still has to solve the problem of Oq = y det || becoming singular for certain values 
of t. To overcome this difficulty one has to glue together various local solutions of the form 
( p.l3| ). Here the crucial observation is that whenever an eigenvalue of || changes sign 

(causing |^ to become singular) this results in a phase jump of the wave function by — |. 
For a closed curve C the number of times this happens along C is a topological invariant of 
C, its Maslov index fj,, see [0, |T6|. On a d-torus T'^ let us choose a set of basis cycles {Cj}, 
j = 1, . . . ,d, such that along Cj the angle -dj increases by 27i. and all other angles ^k^^j 
remain constant. See figure [l| for an illustration of this basis for a 2-torus. Then every 
closed curve on T'^ is a linear combination of the basis cycles Cj. 

All we have to do now in order to get a globally well-defined WKB wave function is to 
make sure that iPwkb{x) returns to its initial value when we follow its value along Cj. In 
other words the phase change along a loop has to be an integer multiple of 27r, i.e. 

If 71 

PdX - i^j- = 2TTnj , rijEZ. (2.22) 



hJc 2 

These are the quantisation conditions ( |1.5|) . Since the action variables / are given by 

Ij = ^ (( PdX (2.23) 



2tt 
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Figure 1: Sketch of a basis {Ci,C2} of loops on a two-torus T^. 



the EBK energies read 

Er, = H(^h(n + ^')^ (2.24) 

where H is the Hamiltonian transformed to action and angle variables and fij denotes the 
Maslov index of the basis cycle Cj. Further restrictions on the values which the integers n 
may assume usually arise from the possible values the action variables / can take in the 
particular problem of classical mechanics, cf. the examples given below. 



3 The semiclassical wave function for the Dirac equation 

In this section we review the determination of semiclassical wave functions for the Dirac 
equation and introduce the relevant notation for the following sections. 

The first steps towards a semiclassical wave function for the Dirac equation were un- 
dertaken by Pauli ||20|- He inserted an ansatz similar to ([275| ) into the Dirac equation and 
found that the phase function S has to solve a relativistic Hamilton-Jacobi equation. He 
then solved the analogue of the transport equation ( p.lO| ) for a particular case but did not 
derive a general expression like ( |2.12| ). The problem was taken up again many years later 



by Rubinow and Keller |21| who proceeded one step further. They showed that the solution 



of the transport equation can be related to the Thomas precession [^,[3^, see also of a 
classical spin vector. However, also in this work no general quantisation conditions similar 
to (|1.5|) were given. As we will see in the following sections their construction is com- 
plicated by the occurrence of non-Abelian Berry phases for which additional integrability 
conditions are needed. 

We also list some related literature which we will, however, not directly refer to in 
the following: Semiclassical approximations to the radial Dirac equation were studied in 
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P5| , |36| , |37| . The semiclassical time evolution of the Dirac equation was examined in p8[ . 
Semiclassical quantisation of subspectra of the Dirac Hamiltonian based on the complex 
germ method was discussed in p9l |40|. The time evolution of semiclassical Wigner functions 



for the Dirac equation is addressed in |^ 

We briefly repeat the basic steps in the derivation of semiclassical wave functions for 
the Dirac equation. For details we refer the reader to ||2l| or to ||4^ where the notation is 
similar to that used here. 

The aim is to find asymptotic solutions to the stationary Dirac equation 



with Dirac Hamiltonian 



D 



ca ( yV - ^M^) ) + P'mc^ + 



(3.1) 



(3.2) 



in the semiclassical limit h 0. The wave function is now a four-spinor, \E' G L^(M^) ® C, 
and the 4x4 matrices cx and /3 are given by 



a. 



cr 

(T 



/3 



l2 







-12 



(3.3) 



where each entry is to be understood as a 2 x 2 matrix. The Pauli matrices cr are given 
by 



1 

1 



0-2 



-i 

1 



0^3 



1 
-1 



(3.4) 



and 1„ denotes the n x n unit matrix. The Dirac equation describes a particle of mass 
m and charge e moving under the influence of the external electro-magnetic potentials 
{(t>{x), A{x)), i.e. we have fixed a frame of reference in which the potentials are static. 

We now modify the semiclassical ansatz (|2.5| ) such that the amplitudes ak{x) take values 
in but the phase S{x) is kept scalar. Inserting this ansatz into the Dirac equation ( p.l| ) 
yields 



[Hb{V^S, x)-E]ao + j{ [Hr){V^S, x) -E]ai + ca(V.ao) } + (^') = 

where Hd{p,x) denotes the Weyl symbol of the Dirac Hamiltonian (|3^), 

mc^ + e(j){x) cr {cp — eA{x))^ 



(3.5) 



Hd{P,x) 



mc — e0[x) 



(3.6) 



cr (cp — eA{x)) 

When comparing ( p.5| ) with ( |2.6| ) notice that ca. = ^ifo- Since the leading order equation 

[HT){V^S,x)-E]ao = (3.7) 
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is a matrix equation it implies the necessary condition that the expression in square brackets 
has an eigenvalue which vanishes identically, i.e. 

H^{V^S,x) = E, (3.8) 

where the eigenvalues 

H^{p, x) = e(l){x) ± ^{cp- eA{x)y + m'^c^ (3.9) 

of Hd act as classical Hamiltonians for our problem. In the following we will use the 
notation S"^ indicating which of the two Hamilton- Jacobi equations (|3^) is solved by 



the respective phase function; the solution itself is again given by (|2.9| ) where now the 
integration is along the flow lines of . Each of the eigenvalues j ^M) has a multiplicity 
of two and the corresponding eigenvectors can be chosen as the columns of the 4x2 
matrices V±, 



V,{p,x)= ' f e + mc^ 



>y2e{e + mc2) \(^icp - eA{x)) 
, . 1 f(7{cp- eA{x)) 

V-ip, X) = . , , 2\ 



(3.10) 



where we have introduced the abbreviation 

e{p, x) := ^{cp - eA{x)f + m^c^ . (3.11) 

With this choice the eigenvectors are orthonormal and complete, i.e. 

Vl{p, x) V±{p, x) = l2, V^ip, x) r±(p, x) = 0, 
and V+{p, x)Vl{p, x) + V^{p, x)Vl{p, x) = U. 



(3.12) 



Now 5"^ (a;) solving one of the Hamilton- Jacobi equations (|3.8| ) is not sufficient in order to 
satisfy eq. ( ^3.71) . Instead we also need ag (a;) to be of the form 

a^{x) = V±{V^S,x)b^{x) (3.13) 

with the still unknown function b^{x) taking values in C^. 

An equation for can be obtained from the next-to-leading order equation deriving 
from ( |3.5|) by inserting ( p.l3| ) and multiplying with V^(VxS,x) from the left yielding 



cViiV^S, x)cxV{VxS, x) (V,6±) + cVliV^S, x)ct[{VxV){VxS, x)] b^ = . (3.14) 
After some algebra one arrives at 

[{VpH^){VxS,x)] (V,6±) + ^ [V,(Vpi/±)(V,^±,x)] 6± + '-aB^{\/xS^,x)b^ = 

(3.15) 

with B^(p, x) = T—B(x) + — [cp - eA(x)] x E(x) . 

e e{e + mc') 
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Here we have introduced the electric and magnetic fields 

E{x) = -V(f){x) and B{x) = V x A{x) . (3.16) 

Equation ( |3.15|) is of the same form as the scalar transport equation ( |2.10| ) except for the 
extra term ^crB^ h^. Thus we can exploit our knowledge of how to solve (|2.10|) by making 
the ansatz 

h^{x) = ^det^u^ix) (3.17) 
leaving us with the spin transport equation 

ii^ + ^aB^{Va,S^,x)u^ = (3.18) 

for the C^-valued function u'^lx), where the dot denotes a derivative along the Hamiltonian 
flow 0*j:^±, cf. eq. ([21I| ). 

Given an initial value u^{y) at a point y its value at a point x, which is connected to 
y by the trajectory (£,?/), is given by 

u^{x)=d±{ty,t)u^{y), (3.19) 

where (i±(^, t/, t) is a 2 x 2 matrix. We have explicitly indicated the dependence of d± on 
the initial point in phase space (^, y) where we start the integration and the time t until 
which we proceed. Clearly, d± also has to solve a spin transport equation, 

d±{^, 2/, t) + \(tB^{<P'h± y)) d±{t y,t) = 0, rf±(C, y, 0) = I2 . (3.20) 

Since the coefficient |crB^ takes values in the Lie algebra su(2) it follows that d±{^, y, 0) G 
SU(2). More precisely, the matrix valued function d± is an SU(2)-valued cocycle of the 
flow as one easily verifies the composition law 

d±{t y, t + t') = d±{ct>'H±{t y),t') d±it y, t) . (3.21) 



Accordingly we can define the skew product fiow |^ 



Yl: M'^ X R'^ X SU(2) ^ M'^ x M"^ x SU(2) 

, , , (3.22) 

{P,x,g) ^ [(p^±{p,x),d±{p,x,t)g) 

which preserves the product of Liouville measure on phase space M°' x M"' and Haar mea- 
sure on SU(2), see |^. Since the cocycle d± takes values in the group G = SU(2) this 
construction is also known as a group extension or, more precisely as an SU(2)-extension. 

At this point the spin degree of freedom is still described on a quantum mechanical 
level in the sense that it is represented by elements of which are evolved in time by an 
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SU(2)-valued propagator. It is, however, possible to switch to a purely classical description. 
To this end consider the adjoint representation of SU(2) defined by 



Ad 



9 ■■ 5u(2) 
Z 



su(2) 

9Z9-' 



(3.23) 



Since an element of the Lie algebra su(2) can be written as a linear combination of the 
Pauli matrices, Z = za, z G M.^, expanding also the right hand side of (|3.23|) in this basis, 
gZg~^ = (99((yf)z)cr, provides us with a rotation matrix (p{g). The map ip : SU(2) — > S0(3) 
is two-to-one and known as the covering map. Anticipating classical spin as a vector of 



constant length let us consider sq 
solves 



On easily verifies that s = ip{d{^,y,t)) Sq 



^ = B±(0*^4^,y))x5 



(3.24) 



i.e. the equation of Thomas precession 32, p3|, which has thus been derived from the Dirac 



equation 0, |4T|. Classical spin precession ( p.24| ) and the Hamiltonian flow 
can now be combined into the classical skew product ||45| 



yt 
-' cl± 



X 



{p, X, s) 



(3.25) 



As opposed to this is a symplectic flow conserving the product of Liouville measure on 



the phase space 



of the translational degrees of freedom and the surface element 



on the sphere which acts as the phase space for the spin degree of freedom. As we will 
see below the properties of determine the semiclassical solutions of the Dirac equation 
as the properties of determine that of the Schrodinger equation. 

Collecting the results of the previous paragraphs, so far we have obtained the following 
expression for the semiclassical wave function for the Dirac equation. 



det et^*(-) V^iV^S^, x) rf±(^, y, t) u^{y) 



(3.26) 



The new ingredients as compared to ( |2.13| ) are the projection matrices V± selecting either 



positive or negative kinetic energies and the spinor part d±{^,y,t)u (y). We can now 



explain why ( p.26| ) does not immediately lead to a generalisation of the EBK quantisation 
conditions (|2.22|) . Assume that the Hamiltonian flow generated by either H~^{p,x) or 

H~ {p, x) is integrable in the sense of Theorem [l|. Then we know that the factors y^det|| 

and e^'^''''^'") can be deflned globally by the construction described in the preceeding section. 
This also applies to the projection matrix V± which simply has to be evaluated along the 
respective path of integration; in particular, for a closed path also the value of V± returns 
to its initial value thus not contributing to the quantisation conditions. 

In contrast, integrating the spin transport equation (p.20| ) along a closed path Cj will 
yield an SU(2) matrix dj. Thus the initial and flnal values of the semiclassical wave function 
( p.26| ) would not just differ by a phase factor but by an SU(2) transformation which would 
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have to be compensated for in order to turn \l/sc(a^) into a globally single-valued object. 
Even worse, integration along a different loop C^^j (which we still have to define, so far we 
are only able to integrate (p.20|) along the fiow lines of 0^±) will in general yield a different 



SU(2) matrix dk=/=j which need not commute with the first one. This can make it impossible 
to find a globally well-defined semiclassical solution as was pointed out by Emmrich and 
Weinstein in a general setting p^]. From the point of view of physics this is not surprising. 
We know that in order to be able to use EBK quantisation the corresponding classical 
system has to be integrable. Now we are dealing with a system that has an additional 
degree of freedom, namely spin. So far, however, we have not imposed any condition on 
the spin dynamics but we have only required the Hamiltonian flow to be integrable. 
In the following section we will show how the notion of integrability can be extended to 
skew products of the form (|3.22|) or (|3.25|) and prove a generalisation of Theorem 



4 Integrability of skew products 

From the point of view of semiclassics the crucial aspect in the characterisation of integrable 
systems by Theorem |l| is the geometrical description of the invariant manifolds which is 
a consequence of the existence of commuting flows (f)\, . . . , 0^. Therefore, one has to find 
a generalisation of the condition We will show that in this sense the following 

definition provides a good generalisation of the notion of integrability. 

Definition 2. The skew product Y^^^ is called integrable, if 

(i) the underlying Hamiltonian flow (j)^^± is integrable in the sense of Liouville and Arnold 
(Theorem^ , i.e. besides the Hamiltonian H^{p, x) =: Ai{p, x) there are d—1 more 
independent integrals of motion, A2{p,x), . . . , Ad{p, x) with 

{A,,Ak} = \/ j,k = l,...,d (4.1) 

and 

(a) the flows 025---)0d can also be extended to skew products l^i* on M'^ x M"^ x 
fy^ii = ^di) with fields Bj{p,x), i.e 

l^i5(p, X, s) = (0* (p, x), ifidjip, X, t))s) (4.2) 

dj{p, X, t) + ^<tBj((/)*(p, x)) dj{p, x,t)=0, dj{p, x, 0) = I2 , (4.3) 

fulfilling 

{A„Bk} + {B„A,}-B,xBk = ^ j,k = 1, . . . ,d . (4.4) 

In view of what was said above let us first show that this is indeed a good definition by 
the following lemma. 
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Lemma 3. Two skew products X-i* and Y^l. of the type ( fl.2| ) commute if and only if the 
corresponding base flows 0*- and commute and if the fields Sj and 3k fulfil 



{Aj, Bk} + {B„ Ak} -BjxBk = 0. 



(4.5) 



Proof: Since the two skew products can only commute if the corresponding base flows 
commute it remains to show that 



for all (p, X, s) eW^ xR'^ X S'^ and all t, t' G M, or equivalently 

'^{dj{(pi{P,x),t')) ^{dkip,x,t)) = (/?(4(0*'(p,a3),t)) ^{dj{p,x,t')) 



(4.6) 



(4.7) 



for all (p, aj) G M*^ X W'- and all t, t' G M. Moreover, since is a double covering, ( |i75| ) is 
also equivalent to 



dj{(t>liP, x),t') dk{p, X, t) 4((/'*'(p, x),t) dj{p, X, t' 
However, due to 

dji^iip, x),t') dk{p, X, t) dk{(t)^-{p, x),t) dj{p, X, t' 

we conclude that 



±l2. 



t'=0 



^ A{t,t') = 0, 
where the difference A{t,t') is defined by 

A{t,t') := dj{(f)i{p,x),t') dk{p,x,t) -dk{(p*-{p,x),t)dj{p,x,t') . 

It is easy to see that A{t, 0) = A(0, f) = V t, t' and thus 

A(t, t') = ^(0, 0) tt' + (t^ + t'2) , t,t'^0. 
The relevant second derivative is given by 



d^A 



dtdt 



7(0,0) 



d_ 

di 



-^aBj{(f)l{p, x)) dk{p, X, t) 



t=o 



d_ 
dt' 



(tB, 



{p,x))dj{p,x,t') 



J t'=0 



■^(T{Ak,Bj}{p,x) + ^{ctBj{p,x)) U(TBk{p,x)) 



'^(T{Aj,Bk}ip,x) + U(TBk{p,x)) ^{(TBj{p,x)) 



)-a Au} + {Aj, Bk} - Bj x Bk] (p, x) 



(4.8) 

(4.9) 
(4.10) 
(4.11) 

(4.12) 



(4.13) 
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which already proves one half of Lemma |^: If the skew products Ydj and X-i^ commute 
then (|4.5|) holds. For the reverse direction notice that (|4.5|) now implies 

A{t,t') = {t^ + t'^), t,t'^0. (4.14) 

Dividing the time intervals t and t' into N subintervals of length e = t/N and e' = t'/N, 
respectively, we can rewrite the first term in (^.11| ) as follows, 



dj{(l>k{P^x),t') dk{p,x,t) 

= dj{(f)l o x),t' - e') dj{(pl{p, x),e') 4(0r'(p, x),e) dkip,x,t- e) 



= :(*) 



(4.15) 



Repeating this procedure A^^ times yields 

A{t,t') = N^e^ + e'^)=N^l/N^) = {l), N ^ oo , (4.16) 

i.e. A(t,t') vanishes which proves Lemma ^. □ 
Having thus established a reasonable generalisation of the notion of integrability we 
can now state the main result of this section. 

Theorem 4. // the skew product flow Y^y± integrable, the combined phase space R'^ x 
X S"^ can be decomposed into invariant bundles Tq T'^ over Liouville- Arnold tori 
T'^ with typical fibre . The bundles can be embedded in T'^ x S"^ such that the fibres are 
characterised by the latitude with respect to a local direction n{p,x), i.e. 

Te = {(p, X, s)eT''xS^\ <(s, n(p, x)) = 6} . (4.17) 

As we have seen above, integrability of Y^-^ also implies similar properties of F*. Let 
us therefore state the following proposition for group extensions, the proof of which will 
facilitate the proof of Theorem ^. 

Proposition 5. Let be an integrable Hamiltonian flow in the sense of Theorem A 
group extension Y^ of (f)\ with group G = \J{n) or a subgroup thereof, 



Yl : X M.'^ X G ^ xM.'^ xG 

{P,x,g) ^ {(p\{p,x),d{p,x,t) g) 



(4.18) 



d{p,x,t) + M{(f)\{p,x)) d{p,x,t) = 0, d{p,x,0) = ln, M : R'^ x R"^ ^ g , (4.19) 

where g is the Lie-albegra of G, is called integrable if the flows can also be 

extended to G-extensions Y2, ■ ■ ■ ,Yj with 

{A,,Mj + {M„A.} + [M,,Mfc] = Vj,fc = l,...,rf. (4.20) 

Then ( ^1.19| ) defines a connection in the trivial principal bundle T'^ x G whose holonomy 
group is an Abelian subgroup of G. 
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We remark that in the case of G = SU(2) the matrices Mj take the form 



Mj{p,x) = -(TBj{p,x) 



(4.21) 



thus reducing condition ( |4.20| ) to 
Proof of Proposition ^: First notice that condition (|4.20|) ensures the commutativity 
of the G-extensions Y^, . . . ,Y^. This can be seen by repeating the proof of Lemma ^ where 
in the definition (|4.11| ) of the difference A{t,t') the matrices dj and dk take values in G, 
see ( [1.19| ). Then equation (|4.12|) is still valid and ( [1.13D reads 



dtdt 



;(0, 0) = {{Mj, Ak} + {A„ Mk} + [M„ Mk]) (p, x) . 



(4.22) 



The following steps up to can be adopted identically. 

Now define the multi-time fiow 

Y* := Y^'o-.-o y/i , 



(4.23) 



where due to commutativity of the fiows Y^^ , . . . , Y^'' ordering is irrelevant. Explicitly we 
have 



Y\p,x,g) = {{(f)'^' o...o(j){^){p,x), d{p,x,t)g) with 
d(p, x,t) : = dd (^{(Pd-i ° ■ ■ • ° 4>i){P, x),td^ ■ ■ ■ di{p, x, ti)g . 



(4.24) 
(4.25) 



Consider a Liouville- Arnold torus T"^ which is invariant under the restriction of Y* to 
M.'^ X M''. The multi-time cocycle ( [1.25| ) then defines a connection in T'^ x G. In order to 
determine the holonomy of this connection choose a basis {Cj} of closed loops on T'' as 
described in section ^ With each loop Cj we can associate a unique (minimal) tuple tj 
such that 

$*(p,a;), tfc G [0, (t,)fc] , k = l,...,d (4.26) 

topologically describes Cj for any (p, x) G T'^. We denote the cocycle (d(p, x, tj) associated 
with Cj by dj{p,x), i.e. 

Y'^{p,x,g) = {p,x,dj{p,x)g). (4.27) 

We immediately see that due to the commutativity of the skew products Yj the cocycles 
dj{p,x) commute thus generating the Abelian subgroup 



H, 



geG 



g = ll[dj{p,x)p ,neZ' 



(4.28) 



of G. In order to see how two such groups, say H{^p^x) and H(^p>^x'), are related, recall that 
for any two points {p, x) and {p', x') on a Liouville-Arnold torus there exists a d-tup\e t 
such that 

^\p,x) = ip',x'). (4.29) 
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Again due to commutativity we have the equahty 



Y*^- = Y"* o Y*^- o Y* (4.30) 

dj{p, x) = d(p', x', -t) dj{p', x') d(p', x', t) . (4.31) 
Y-*oY* = id => d{p',x',-t) = [d{p',x',t)]-^ (4.32) 

H{p',a^') = gH[^p^^)g-^ , g := d(p', x', t) , (4.33) 

i.e. the subgroups H(^p/ x') at different points are obtained by conjugation with a group 
element g G G. Thus H{^ x) and H(^pi x') are isomorphic and we have identified the Abehan 
group ( [4.28|) as the holonomy group of the connection defined by ( |4.19| ). □ 



implying 
Moreover, 
and thus 



Proof of Theorem |: Applying Proposition || to the SU(2)-extension (|3.22|) , with each 



point (p, x) on a Liouville-Arnold torus is associated an Abelian subgroup H(^p x) of SU(2). 
Abelian subgroups of SU(2) are either one-parameter subgroups or discrete subgroups 
thereof, i.e. we can associate with each point (p, a:;) G T"^ a one-parameter subgroup of 
SU(2). The latter can be parametrised as 

Hr, := {g G SU(2) | g = e-^t"^" , a G [0, 47r) } (4.34) 

with a direction characterised by the unit vector n & S'^ . By means of the covering 

map if, see definition below (|3.23| ), this construction uniquely determines a one-parameter 
subgroup of S0(3), ^{Hn{p,x)) , at each point (p, cc) of a Liouville-Arnold torus. This fact 
in turn allows for a construction of invariant manifolds of F^*. Consider a point (p, x) G T*^ 
and a spin vector s E S"^ ^ M?. Transporting s along a path on T"^ by means of the 
multi-time fiow 

Y*i := Yj^^o...o Yj,^ : (p, x, s) ^ ((</>*'' o ■ ■ ■ o </,*/)(p, x), vp(d(p, x, t)) s) (4.35) 

gives rise to rotations (y9(d(p, cc, t)) of s. The rotation associated with a closed path Cj is 
given by the rotation matrix ip{dj{p,x)) G ip{Hn{p^x)) C S0(3). Thus whenever the path 
on the Liouville-Arnold torus is closed the spin vector s returns to a point on the circle 
ip{Hn[p^x))s, which is a parallel of latitude with respect to the axis n(p, x). Corresponding 
circles at different points (p', x') are obtained as follows, see figure |. The one-parameter 
subgroups at different points of the torus are related by conjugation with d(p, a;,t), see 
( [1.31J ). However, from the definition of the covering map ip, see below (|3.23|) , it follows 
that 

Hn{p',x') = <d(p, X, t) Hn{p,^) [d(p, X, t)]"^ = -f/'^(<ri(p,cc,t))ri • (4.36) 

On the other hand from 

Y*i(p, X, t) = (p', x\ ^(d(p, X, t))s) (4.37) 
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Figure 2: Illustration of the invariant manifolds Tq, see ( 4.17] ). At different points of the Liouville- 
Arnold torus T*^ spin vectors are restricted to parallels of latitude with different axes n; the angle 
<{s,n{p,x)) = <{s' , n{p' , x')) is conserved 



we see that by moving from (p, x) to (p', x') the spin vector s is rotated in the same way 
as is the axis n. Therefore the angle 

9 := <{s,n{p,x)) (4.38) 

is conserved by Y*i and thus by Y^i, and the bundles Tg, see (|4.17| ), are invariant, concluding 
the proof of Theorem 0. □ 



5 Quantisation and spin rotation angles 

With the novel notion of integrability for skew products at hand we can now return to the 
semiclassical wave function ( |3.26| ) of the Dirac equation. We have already seen that except 
for the spin transporter d±{^,y,t) all terms in the local expression ( |3.26| ) can be given 
a global meaning provided that the classical translational dynamics is integrable. In the 
case of the Schodinger equation, this observation led to the EBK quantisation conditions 
in a straightforward way, see section |[ We will now show that this is also the case for the 
Dirac equation if the skew product flow 1^*^ is integrable in the sense of Definition ||. 

Theorem |J provides us with a method for integrating the spin transport equation not 
only along the flow lines of but also along those of the flows (pl, . . . , 0^, which are 
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defined in the Theorem [l|. The resulting spin transporter is defined on the whole Liouville- 
Arnold torus and given by the multi-time cocycle <d{p,x,t), see (|4.25|) , for a point {p,x) 
which is reached from {$,, y) by $*(^, y) = {p, x). For a closed path Cj the initial and final 
values of the spin part of the semiclassical wave function are still related by the SU(2) 
transformation dj{$,,y). However, provided that the skew product 1^*^ is integrable we 
know that the matrices dj{$,,y) for different loops Cj commute. Therefore we can choose 

to be a simultaneous eigenvector of these SU(2)-matrices which effectively reduces the 
SU(2)-holonomy to a simple phase, i.e. a U(l)-holonomy. 

Introducing the local axis n(^, y), see Theorem ||, we have 

^i(^,2/)=e-^"^'^"(«'^), (5.1) 

where aj is the angle by which a classical spin vector is rotated about the axis n(^, y) when 
transported along Cj; note that aj has to measured modulo 471. The spin part u^{y) of 
the wave function ( |3.26|) is chosen to be an eigenvector of crn(^, y) with eigenvalue either 
+1 or —1. Thus the phase shift resulting from spin transport is given by e^^""]/"^ . Together 
with the contribution from the classical action and the Maslov phase, cf. ( |2.22|) , the total 
phase shift suffered by the wave function when going through the cycle Cj is 

i£pdX-^,,Tf. (5.2) 

Requiring this to be be an integer multiple of 2-n in order to obtain a single valued wave 
function yields the novel quantisation conditions 

(b PdX = 27ih (uj + !^ + ms—] (5.3) 
Jcj ^ 4 271/ 

where we have introduced the spin quantum number = ±|. Equation (|5.3|) is our 
central result replacing the EBK quantisation condition ( |1.5| ) in the case of relativistic 
particles with spin |. If the classical Hamiltonian H^{p,x) expressed in action and angle 
variables see ( |2.23|) , is given by H{I) then the semiclassical eigenvalues resulting 



E^^^^=H(h(n + ^ + m,— )) . (5.4) 



from (|5.3|) read 

4 271. 

So far we have only employed the fiows Y^j resulting from the "Hamiltonians" Ai = 
if ^, A2, . . . ,Ad complemented with fields Bj{p, x). For practical purposes it is helpful to 
directly make use of the fiows generated by the action variables / in the semiclassical 
quantisation process. To this end these fiows have to be extended to skew products on 
M"^ X M'^ X 5^ by some suitable fields Bi. . Transport along a basis cycle Cj is then given by 
Y^J" and the relevant SU(2) transformation dj = exp{—iaj(Tn/2) is given by the cocycle of 
just one fiow instead of the a linear combination of the d fiows. Since 



20 



where cj are the fundamental frequencies, see ( p. 17] ), we have to require that 



(5.6) 



yielding the consistency condition 



(5.7) 



We will illustrate these remarks when applying the method to explicit examples in sections 
0-1- 



6 Non-relativistic limit: The Pauli equation 

In this section we show that the semiclassical quantisation scheme developed in the pre- 
ceding sections for the Dirac equation carries over to the Pauli equation. The latter arises 
as a non-relativistic approximation to the Dirac equation but can also be generalised to 
describe particles with spin s G No/2 other than s = |. We show that also in the case 
with s 7^ I the semiclassical analysis of the Pauli equation gives rise to a skew product 
on M*^ X M*^ X S"^ which can then be quantised along the same lines as in the case of the 
Dirac equation. We keep the presentation short relying heavily on our treatment of the 
Dirac equation in the preceding sections; details on the semiclassical analysis of the Pauli 



equation with arbitrary spin can be found in |^ 

The Pauli equation for a spin-| particle can be obtained from the Dirac equation ( pT 



in the non-relativistic limit c — > cxd, see e.g. |4^. With the representation ( 3.3|) of the 
matrices ot and /3 and writing the Dirac spinor as = {ip'^, x^) with ip,x ^ L^(]R^) (g) 
one finds the following equation for the upper two components, 

Hp^{x) = Eij{x) (6.1) 



with Pauli Hamiltonian 

Hp = -i—{-V--A{x)\ +e<j){x)-—B{x)'-(T. (6.2) 
2m \ 1 c / mc 2 

Let us generalise our discussion to Pauli Hamiltonians 

^P = i^s + ^d7r,(o-)B, (6.3) 

where Hs is a Schrodinger Hamiltonian with Weyl symbol H{p, x), dusicr) denotes the 2s + 
1 dimensional (derived) irreducible representation of 5u(2) and B is the Weyl quantisation 
of the classical vector valued function B{p, x) on phase space. The wave function has 
now 2s + 1 components, i.e ip ^ L'^iM.'^) ® C^'^"*'^. The special case (|6.2|) is recovered by 
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the choices d = 3, s = ^, H{p, x) = ^ {p - -^A{x)Y + e(f){x) and B{p, x) = -^B(x). 
Spin-orbit coupling can, e.g., be described by 

BUp,x) = f{r)L, (6.4) 

where L = x x p is orbital angular momentum and / is an arbitrary function of the radial 
coordinate r = \x\. 

As in (|275|) we make the ansatz 

^(^) = E(t) ^k{x) e'^'^^'^ (6.5) 

with scalar S, and taking values in C^*+^. Upon inserting ( |6.5| ) into the Pauli equation 
( |6.1|) with Hamiltonian (|6.3| ), in leading order we find the Hamilton- Jacobi equation 

H{V^S,x)=E (6.6) 

and in next-to-leading order we obtain the transport equation 

[{VpH){V^S, x)] (V.ao) + ^ [V,(Vpi/)(V,5, x)] + ^d7r,(o-)B(V,5', x) = Q . (6.7) 

As below (|3.15|) the ansatz 



dy 

aQ{x) = \l dei — u{x) (6.8) 
V ox 

leaves us with the spin transport equation 

u{x) + ^dvr, (o-)B( Vcc^, x)u{x) = {). (6.9) 

Integration along a trajectory with 4>^ii{^,y) = {p,x) yields 

u{x) = 7cMty,t))u{y), (6.10) 

where vr^ denotes the 2s + 1 dimensional unitary irreducible representation of SU(2) and d 
solves the spin transport equation 

d{t 2/, t) + ^o-e(0*,(C, y)) dit y,t) = 0, d{^, y, 0) = l2 , (6.11) 

in SU(2). Thus the covering map ip : SU(2) SO (3) still relates spin transport to classical 
spin precession, 

s = B{<P'H{ty))xs, seS^^R\ (6.12) 
and we identify the skew product flow 

Y^,{p,x,s) = {^'HiP,x),MP,x,t))s) (6.13) 
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as the classical system corresponding to the Pauli equation with Hamiltonian (|6.3|). If Y^j 
is integrable in the sense of Definition H we can use the same construction as in section H to 
define a semiclassical wave function associated with a Liouville- Arnold torus T"^. Locally 
we have 

ij^x) ~ \/det||e^^(^) TiM^, y, t)) u{y) . (6.14) 

and the initial and final values of u after transport along along a basis cycle Cj are related 
by TTsidj^^jy)). As in ( |5.1| ) we obtain 

rf,(^,t/)=e-i"^'^"(^'^), (6.15) 

where the axis n is specified by Theorem ^ and the angle aj is the rotation angle for a 
classical spin vector. The representation matrix 7is{dj{^,y)) has eigenvalues exp(— im^aj), 
rris = —s,—s + l,...,s, see e.g. Upon choosing n{$,,y) to be an eigenvector of 

Tisidji^, y)) the total phase change of experienced by -^sc along a cycle Cj is given by 

1 / TT 

— a) P dX — i—^j — nisOij (6.16) 
" Jcj 2 

thus resulting in the semiclassical quantisation conditions 

PAX = 2nh (uj + ^ + rris^^ , (6.17) 

with n E Z'^ and = —s, —s + l,...,s. 



7 Spherically symmetric systems 

Before we treat some explicit examples in order to illustrate the novel quantisation con- 
ditions (|5.3|) and ( |6.17|) we show how to apply them to an important class of integrable 
systems, namely spherically symmetric systems. 

A spherically symmetric Dirac Hamiltonian has the structure 

h 

Hb = ca + p mc^ + e^{r) , (7.1) 

where the electrostatic potential depends only on the radial variable r := |a:;|. One 
easily verifies that the Hamiltonian ( [7.1|) commutes with all components of total angular 
momentum, 

Since one also has 

[Jd, Jbz] = (7.3) 



23 



- 2 

one can choose the eigenfunctions of Hj) to be also simultaneous eigenfunctions of and 

Jdz- 

Analogously a spherically symmetric Pauli Hamiltonian is of the form 

= -|^- A + e0(a;) + hn^icr) f{r)L , (7.4) 
cf. ( |6.4|) , and commutes with all components of total angular momentum 

Jp :=L + ^d7r,(o-). (7.5) 

Again due to 

[Jp, Jpz] = (7.6) 
the eigenfunctions of ^p can be chosen such that they are simultaneous eigenfunctions to 

- 2 

the modulus squared Jp and the 2;-component Jpz of total angular momentum. 

In the semiclassical limit the Hamiltonians (|7.1J ) and ( |7.4D both give rise to skew prod- 
ucts on R'^ X M'^ X S*^ with a spherically symmetric classical Hamiltonian H{p,r) and 
spin precession with fields of the form 

B{p,x) = f{r)L. (7.7) 

It is well-known that 0^ is integrable in the sense of Liouville and Arnold since 

{H, L} = {H, M} = {L, M} = (7.8) 

where L and M are the modulus and 2;-component of classical orbital angular momentum 
L = X X p, respectively. We remark that L{p,x) not being smooth at points where p\\x 
will not play a role in the follwing since all relevant constructions will stay away from these 
points. Thus (|7.8|) holds wherever needed. 

In order to show that Y^^ with field ([7.7|) is also integrable we have to extend the fiows 
0^ and (t)\j to skew products with fields 13 l and 3m fulfilling ([4.4|) . To this end consider 
the Weyl symbol of the operator J := L + |cr, 

J = pxx + ^cr. (7.9) 
By straightforward calculation one finds 

J ■=\j\ = L + -aY + {h'^) and J, = M + -a^ . (7.10) 



Comparing these to the Pauli Hamiltonian (|7.4| ) suggests that 



Bl:=j and Bm := e, , (7.11) 
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where is the unit vector in z-direction, might be a good choice. As one easily checks the 
skew products 1^*, ¥[ and obtained in this way indeed form a set of commuting flows 
on M'^ X M*^ X S"^ and thus F/j is integrable, cf. Theorem |[ 

Spherically symmetric Hamiltonians H{p,r) can be separated in spherical coordinates 
{r,9,(j)), see e.g. ||31|. Introducing the action variables 

Ir = ^ Prdr , Iq = j) p0 and = ^ p^dcj) , (7.12) 

where Pr, Pe and p^p denote the canonical momenta conjugate to r, 6 and 0, respectively, 
one always finds 

I^ = M and h = L - M . (7.13) 

The radial action depends on the particular system under investigation and solving 
Ir = § Pr{E, Iq, I^)dr for E yields the Hamiltonian H in action and angle variables. It is 
a general feature of spherically symmetric systems that H does not depend on and 
independently but only on the sum L = Iq + I^, i.e. if is a function of Ir and L. 

The motion generated by M = J,^ is a rotation about the z-axis and thus there is no 
turning point in the time evolution of the coordinate yielding the Maslov index 

/^M = /i</, = 0. (7.14) 

On the other hand the motion generated by Ig takes place between two turning points for 
the coordinate 9 and thus 

= 2 and accordingly /x^ = 2 . (7.15) 
In order to find the spin rotation angle om one has to integrate the spin precession equation 

s = e^x s (7.16) 

over one cycle of The latter simply changes the angle variable conjugate to M 
with unit speed and we have 

ttM = Svrle^l = 27r. (7.17) 

Similarly 0^ changes with unit speed and since Bl = L/L is constant along such a 
cycle we also find 

aL = 27r\L/L\ = 27r. (7.18) 

Therefore, for any spherically symmetric Dirac or Pauli equation we have semiclassical 
quantisation conditions 

L = h (l + 1- + rris] and M = h{mi + ms) (7.19) 



2 



with integers I and mi. 



25 



Classical mechanics imposes the restrictions 



L > and M <L. 



(7.20) 



Defining the new quantum numbers j := / + and rrij := mi + nis after some trivial 
algebra one finds that ( |7.20|) translates to 



j > and rrij 



(7.21) 



Notice that if s is integer or half-integer, respectively, then so are j and rrij. Finally the 
semiclassical quantisation conditions for angular momentum read 



and M = hrrij . 



(7.22) 



As mentioned at the end of section ^ it will now be useful if we also find a field Br which 
turns the flow 0* generated by the action variable into a skew product that commutes 
with ¥[ and Ylj. To this end we may exploit the relation ( ^.7|) . Keeping in mind that H 
is a function of Ir and L = Ig + only, i.e. ug = = ul = dH/dL, this yields 



Br 



B - ujlBl 



IjJr 



(7.23) 



Since B = f{r)L is parallel to Bl = L/L and since L is a constant of motion Br does not 
change its direction along the flow line of either 0^ or 0*. Therefore 



ar 



£5 — ujly 

i^dl^r 

UJr 



(7.24) 



where 'dr is the angle variable conjugate to Ir. Splitting the integrand into two terms, 
B/oJr and —ujlL/ (urL)^ the second integration is again trivial. In the first integral we can 
make use of &drlujr = &dr/'dr = dt, where t is the physical time along a fiow line of 05?. 
Therefore we have obtained the handy expression 



ar 



B(0*,(p,a;))dt-27r^^^ 

UJr J-' 



radial 



(7.25) 



where the remaining integral extends over one cycle of the radial motion, e.g. from peri- 
helion to aphelion and back. The missing quantisation condition then reads 

Ir = 2nh (nr + l + rris^] , (7.26) 
y 2 Ztc J 

where we have inserted Hr = 2 for a typical radial motion between perihelion and aphelion. 
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8 Example 1: 

Harmonic oscillator with spin-orbit coupling 

As a first example we discuss the three dimensional isotropic harmonic oscillator with 
spin-orbit coupling in a non-relativistic context with arbitrary spin. 
The Pauli Hamiltonian 

i^P = -|^A + |c.V + ^dvr.(^) kL (8.1) 

describes an oscillator with frequency u and for the spin-orbit coupling we have chosen a 
Thomas term, i.e. the leading order expression of , see ( p.l5| ), in the non-relativistic 
limit c oo with 

eE = -V _ ,2. 

i.e. 



eE = -V y-^^ ) = -muj^x (8.2) 



By comparison with ( |8.1| ) we have k = tu^/ (2mc^) but we may keep k arbitrary for the rest 
of the section. 

One usually determines the eigenvalues of ( |8.1j ) by first observing that besides Jp and 

Jpz, see (|7.5|) , also L and = [h/ 2(171 s{o-)Y commute with Hp and that the spin-orbit 
term can be expressed in terms of these conserved quantities, see e.g. ||50| . 

Obviously, the Hamiltonian ( p.l| ) defines a spherically symmetric system as discussed in 
the preceeding section. Therefore, we immediately have the quantisation conditions ( |7.22|) 



for total angular momentum. The semiclassical eigenvalues, however, can even be found in 
a more straightforward way than to draw on (|7.25|) . To this end notice that the classical 
Hamiltonian H{p,r) = ^ + ycuV^ transformed to action and angle variables reads 

H{Ir,L) =uj{2Ir + L) . (8.4) 

Semiclassical quantisation with s = yields 

En^i = hw (^2n + l + ^^ , (8.5) 

exhibiting the zero point energy of 3hw/2 of the three independent oscillators in x-, y- 
and z-direction. Since H only depends on a linear combination of and L the problem 
has an even higher degeneracy than general spherically symmetric systems. Thus we may 
introduce the new action variable 

/i := 2Ir + L (8.6) 



with corresponding frequency lui = dH/dli = lu. From (|8.5|) we see that without spin Ji 
has to be quantised as Ji = h{2n + / + |). In order to find the correction from the spin 
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contribution we have to extend to a skew product Y^* that commutes with F^*. The 
consistency condition (|5.7| ) uniquely determines the relevant field Si, since it reduces to 

B = uiBi =^ Bi = — = -L. (8.7) 

Integration of the spin precession equation s = Bi x s is once more trivial since X is a 
constant of motion and we get 

ai = \L\ = / + - + mJ. 8.8 

LO LO \ 1 J 

Finally, the semiclassical energies read 

= ^ (2n + / + ^ j + m,f?K n + ^ + m J . (8.9) 



For spin s = ^ these are the exact eigenvalues of the Hamiltonian ( |8.1|) whereas for s 7^ g 
they are good approximations to the exact eigenvalues if / is large, i.e. if the action variable 
L is large compared to fi as required for semiclassical approximations. 

At this point a short remark concerning the practical application of conditions ( ^.31) or 
( |6.17| ) is in order: In general the spin rotation angels a can depend on the action variables 
/. If the dependence is simple as in this example (ai depends on L but not on Ji, and 
au are constant) one can recursively apply ( |6.17| ) as was done above: First we quantised 
L and then used the result when quantising l\. In general, however, one may first have 
to solve the quantisation conditions (|6.17|) for the action variables I before inserting the 
result into the Hamiltonian Hit). 



9 Example 2: 

Sommerfeld's theory of fine structure revisited 

In this section we address the classic problem of Sommerfeld's theory of fine structure 
which was already mentioned in the introduction as a motivating example. 

To this end we have to consider the relativistic Kepler problem with classical Hamilto- 
nian 

e2 



/7(p,r) = -— + VcV + m^c4- (9.1) 

Solutions of Hamilton's equations of motion corresponding to bound states are given by 
"Rosettenbahnen", ellipses with moving perihelia. Since (p.l|) is spherically symmetric 
angular momentum is conserved and the motion takes place in a plane. Introducing polar 
coordinates in this plane the orbits can be expressed as 




' cos , (9.2) 



c2L2 - 



=A 
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where for later reference we have defined the constants C, A and 7. Detailed information on 
the classical mechanics of the relativistic Kepler problem can, e.g., be found in Sommerfeld's 
original article or in his book |51]. Transforming the Hamiltonian to action and angle 
variables yields 



H{Ir,L) 



mc 



1 + 



Sommerfeld quantised the system using (|1.2|) by demanding 

L = hrir and L = hi 



-1/2 



(9.3) 



(9.4) 



with integers and /. Classical mechanics imposes the restrictions /r > and L > 
yielding n^, / > 0. Moreover, Sommerfeld excluded I = because in this case the electron 
would collide with the nucleus. Thus he found the energy levels 



E. 



Sommerfeld 
rirl 



mc 



1 -1/2 



1 + 



at. 



rir eNo, I eN, 



(9.5) 



where as = e^/ (^c) denotes Sommerfeld's fine structure constant. 

When now quantising the relativistic Kepler problem with spin | using conditions ( ^73| 
first notice that due to spherical symmetry we have, cf. (|7.19|) , 



h{ I + - + ms 



2 



(9.6) 



Since the additional conditions (|7.20|) do not allow the combination / = 0, = — | (which 
would lead to a negative quantum number j) we see that L/h E N, i.e. it assumes the 
same values as in Sommerfeld's prescription (|9^) . 

In order to quantise the radial action variable Ir we also have to calculate the spin 
rotation angle to which end we can use the general formula ( [7.25|) . The field B occurring 
in the equation of spin precession can be obtained from B^, see (|3.15|) , by inserting E = 
— V(— e^/r) yielding 



B 



1 



L. 



(9.7) 



e{e + mc^) 

where e was defined in equation (|3.11|) . Without restriction we may choose L\\ez for our 
calculation leading to 



L = -r' 



B 



e + mc^ r 



-(he. 



(9.8) 
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Using E = H{p,r) = —e^/r + e, see (|9.1| ), this expression further simplifies to 
which is nicely adapted for use in ([7.25| ), 



(9.9) 



Bdt 



radial 



radial 



{E + mc'^)r{(f)) + 



E + mc^ 



radial 



E + mc^ + e^/r(0) 



Changing variables to 77 = 70 and inserting ( p.2|) we get 

E + mc^ 



62 

7 



1 - 



radial 



7 



(9.10) 



(9.11) 



where the calculation of the last integral is tedious but elementary. From ( |9.3| ) we easily 
determine the frequency ratio in the second term of (|7.25|) to ul/^t = 7~^ and since we 
have chosen -Llle, we obtain 



_ _ 27r - — 

7 / 7 



27r. 



This remarkable result turns ( [7.261) into the rather trivial condition 



h \ Ur + - + m. 



(9.12) 



(9.13) 



with integer Ur and = ±|. Moreover, classical mechanics demands /r > and thus 
Ir/h G No, i.e. the values assumed are again identical to those predicted by Sommerfeld's 
condition (|9.4|) . 

Finally we have found the semiclassical energies, 

1 -1/2 



E, 



rirlnis 



mc 



at 



2 



at 



(9.14) 



which are identical to those obtained from the exact solution of the Dirac equation 
with potential — e^/r, see e.g. |50l, as can be most easily checked by introducing the 
quantum numbers j and mj of total angular momentum, see below ( [7.2Q ), and the principal 
quantum number n, associated with / := + L, by / = hn. Although Sommerfeld's 
formula ( |9.5| ) also yields the correct energy levels of the hydrogen atom it predicts the 
wrong multiplicities. This problem is rectified by the present treatment as can be checked 
straightforwardly. Consider, e.g., the ground state which in ( p.5| ) is obtained by = 0, 
/ = 1, and thus is non-degenerate. On the other hand in (|9.14|) the same energy is obtained 
by choosing either rir = —1, I = 0, rris = ^, ov rir = 0, I = 1, rris = Alternatively these 
states can be characterised by = 1, j = |, = ±i yielding a multiplicity of 2 as in the 
exact quantum spectrum. 
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10 Conclusions 



In this article we have derived semiclassical quantisation conditions for the Dirac and Pauli 
equations. We have shown that spin yields a contribution of the same order of magnitude 
as the Maslov correction. The spin contribution is determined by rotation angles for a 
classical spin vector that is transported along orbits of the translational dynamics. 

The crucial step in the derivation of the semiclassical quantisation conditions was the 
generalisation of the notion of integrability to certain skew products and group extensions. 
The relevant integrability conditions enabled us to effectively reduce the non-Abelian Berry 
phases, appearing in the analysis of multi-component wave equations, to a U(l)-holonomy, 
i.e. to an ordinary phase factor. The latter can then be incorporated into the quantisation 
conditions. 

We remark that our treatment generalises to arbitrary multi-component wave equations 
for which the principal symbol of the Hamiltonian has eigenvalues with arbitrary but con- 
stant multiplicity. One then has to deal with G-extensions, where G = U(n) or a subgroup 
thereof, of the ray dynamics (generated by the eigenvalues of the principal symbol). Propo- 
sition 1^ provides us with the relevant integrability condition which allows for an effective 
reduction of the holonomy group G to an Abelian subgroup H. The spin rotation angles 
are then replaced by the eigenphases of some unitary representation matrices of H. 

In section |^ we have applied the novel quantisation conditions to the relativistic Kepler 
problem. We have seen that by a freak of nature all relevant spin rotation angles are given 
by 271 and thus cancel (or add up to an integer) with the Maslov term. It is this coincidence 
due to which Sommerfeld was able to calculate the energy levels of the relativistic hydrogen 
atom including spin-orbit coupling 10 years before the Dirac equation was developed. 
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A Wigner-Weyl calculus 

With a differential operator A one can associate an object on classical phase space, its 
Weyl symbol A{p,x), by 

{A^){x) = -i— / f a(p, etP(--) ^{z) d'zd'p. (A.l) 

If is a multi-component object, e.g. G L^(R°') (g) 0^"*+^, then A{p, x) is matrix valued. 
Reverting this reasoning, one can also associate an operator A with a more general symbol 
A{p,x), which does not necessarily correspond to a differential operator, via ( |A.1| ). This 
procedure is known as Weyl quantisation and certain properties of symbols translate to 
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properties of the operators, thus leading to so-called pseudo-differential operators, see e.g. 
1 5^ for an introduction. 

If an operator A can be represented by an integral kernel KA{x,y), i.e. 

(A.2) 



(A.3) 



one obtains its Weyl symbol from 



Inverting this transformation yields 



KAix,y) 



1 



A p 



x + y 



p. 



If the symbol A{p, x) has an expansion in powers of /i, 

A{j),x) = ^n!'Ak{p,x) 



(A.4) 



(A.5) 



fc>0 



the corresponding operator A is called a semiclassical Weyl operator. The leading and 
subleading terms Aq and Ai in the expansion ([A.5|) are known as the principal symbol and 
the subprincipal symbol, respectively. 

Application of a semiclassical Weyl operator to a rapidly oscillating function 

k 



an[x) 



k=0 



ak[x) . 



(A.6) 



is governed by the following theorem; the corresponding statement in a slightly different 
setting can, e.g., be found in (S^, chapter 4.3]. 

Theorem 6. Applying a semiclassical Weyl operator A with a symbol A{p, x) of the form 
( [A.5| ) to a wave function of type ( [A.6| ) yields in leading orders as h 0, 



{A<if,c)ix) = {Ao{Va,S{x),x) ao{x) + - 



AQ(yxS{x),x) ai{x) 



+ iVpAo)iVa,Six),x)Va,aoix) + -aoix)[V4VpAo){Va,S{x),x)] (A.7) 



+ Ai(ya,S{x),x) aoix) 



Notice that the nesting of brackets in the fourth term on the r.h.s. indicates that the 
gradient with respect to p is taken before we set p = VxS, whereas the gradient with 
respect to x is only taken after doing so, i.e. 



For the proof we refer the reader to ||5^, chapter 4.3]. 



d d 



d'^A 



° (V.5(a;),a;)-^(x). 



dxjdxk 



(A.8) 
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